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Inter-Landau-level transitions in the bilayer graphene in high perpendicular magnetic field at the 
filling-lactor v = have been studied. The next-nearest-neighbor transitions, energy difference 
between dimer and non-dimer sites, and layer asymmetry are included. The influence of Coulomb 
interaction is taken into account. The magnetoplasmon excitations in bilayer graphene at small 
momenta are considered within the Hartree-Fock approximation. The asymmetry in cyclotron 
resonance of clean bilayer graphene is shown to depend on magnetic field. At lower magnetic fields 
the energy splitting in the spectrum is due to electron-hole one-particle asymmetry while at higher 
magnetic fields it is due to Coulomb interaction. For the fully symmetric case with half-filled zero- 
energy levels the energy splitting proportional to the energy of Coulomb interaction is found both 
for bilayer and monolayer graphene. 

PACS numbers: 73.22.Lp, 73.43.Lp, 76.40.+b, 78.66.Tr 
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I. INTRODUCTION 

Recent experimental progress has allowed the fabri- 
cation and study of monolayer and bilayer graphene. 
The electronic band structure of these objects is gap- 
less and has a chirality 1 . The monolayer has Dirac-type 
spectrum with linear dispersion and chirality exhibit- 
ing Berry phase it. The bilayer graphene is the unique 
object which combines the parabolic dispersion law of 
quasiparticles near the zero energy point with their chi- 
rality exhibiting Berry phase 2w. This picture is ob- 
tained with the tight-binding Hamiltonian for electrons 
taking into account only nearest-neighbor transitions; the 
one-electron spectrum is symmetric around zero energy. 
Taking into account next-nearest-neighbor transitions re- 
sults in the asymmetry of electron spectrum around zero- 
energy point 2 . Some experimental data demonstrate 
electron-hole asymmetry in cyclotron resonance spectra 
of clean bilayer graphene^— . 

One-particle Landau levels in the bilayer graphene at 
high magnetic fields have been considered in the works 8,9 
taking into account only nearest-neighbor transitions. In 
magnetic field there is a two-fold degenerate zero-energy 
Landau level incorporating two different orbital states 
with the same energy. Taking into account spin and val- 
ley degeneracies, the zero-energy Landau level in a bilayer 
is eight-fold degenerate. For the bilayer with small asym- 
metry there are four weakly split two-fold degenerate lev- 
els near zero energy. The valley and orbital degeneracies 
are lifted, but the electron-hole symmetry is preserved. 

The near-zero-levels are strongly influenced by 
Coulomb electron-electron interaction. The electron- 
electron interaction is an important problem in the exper- 
imental study of cyclotron resonance in monolaye r 10 ' 11 , 
bilayer^J^ and multilayer— graphene, exhibiting 
some properties of a monolayer and bilayer. Interaction- 
induced shift of the cyclotron resonance as a function 
of the filling-factor in bilayer graphene was studied in 



ReLi2,the symmetry breaking in the zero-energy Landau 
level in bilayer graphene is demonstrated ir>i£. 

The charge-density excitations at small momenta were 
considered theoretically within the Hartree-Fock approx- 
imation for monolayer graphen o 19 ' 20 and for bilayer—. 
In the worksS^r— electromagnetic response in graphene 
was calculated numerically in the RPA approximation for 
wide range of excitation momenta. The spin-flip excita- 
tions and spin- waves in graphene were studied in Rcf. 25 . 
In the work o 26 ' 27 intra-Landau level transitions were con- 
sidered. The many-body corrections obtained within 
the renormalization method, including weak electron- 
hole asymmetry, and the attempts to explain sharp tran- 
sition from square to linear dispersion regime were re- 
ported in Refsi^r— . In the worksi^— Coulomb interac- 
tion was shown to conserve electron-hole symmetry for 
excitations. 

In the present paper the inter-Landau-level transitions 
in the bilayer graphene in high perpendicular magnetic 
field at the filling-factor v = are studied. The nov- 
elty of this work is that the electron- hole asymmetry and 
Coulomb interaction are included into consideration si- 
multaneously. Special attention is given to the difference 
in the cyclotron transition energies for two valleys un- 
der different conditions. First, the one-particle Hamil- 
tonian and Landau levels are considered taking into ac- 
count the electron-hole asymmetry due to next-nearest- 
neighbor transitions and energy difference between the 
dimer and non-dimer sites and the small energy difference 
between the layers due to external potential. Then, the 
influence of Coulomb interaction is included. The charge- 
density excitations (magnetoexcitons) at small momenta 
are considered within the Hartree-Fock approximation in 
the case of clean (neutral) bilayer graphene with filling- 
factor v = 0. The energies of excitations are shown to be 
different in the two valleys, and the origin of this differ- 
ence depends on magnetic field. At lower magnetic fields 
the energy splitting is due to electron-hole one-particle 
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asymmetry while at higher magnetic fields the energy 
splitting in the spectrum is due to Coulomb interaction. 
Next, the results are discussed in connection with exper- 
imental possibility to observe the influence of Coulomb 
interaction, and the comparison with other theoretical 
works is presented. 



II. HAMILTONIAN AND LANDAY LEVELS OF 
BILAYER GRAPHENE 

The bilayer is modelled as two coupled hexagonal lat- 
tices with inequivalent sites (Al, Bl) and (A2, B2) in 
the first and second graphene layers, respectively, ar- 
ranged according to Bernal (A2-B1) stacking. In the 
tight-binding model the energy states of electrons in 
(A1-B2) dimcr in the vicinity of zero-energy point are 
conveniently described by an effective two-component 
Hamiltonian2i&2 that operates in the space of wave func- 
tions ^ = (ipAi, i>B?) in the valley K and of f = 
(jpB2,ipAi) in the valley K' . The asymmetry between 
on-site energies in the two layers U arising from the in- 
fluence of external gates or a doping effect, the next- 
nearest-neighbor transitions and the difference between 
on-site energies of dimer and non-dimer sites A are taken 
into account: 



non-dimer sites. The first term is responsible for the 
electron-hole asymmetry in the spectrum around zero- 
energy point. The second term due to A results in the 
shift of the single-energy spectrum as a whole. This term 
is unimportant for our considerations an will be omitted 
later. 

The two-component Hamiltonian is applicable if the 
considered electron energy e\ is within the energy range of 
l £ l < z7i- The weak asymmetry means that U/ji <C 1, 
A/71 < 1, 74/70 < 1. 

In the perpendicular magnetic field B the energy spec- 
trum of Landau levels E ni [n = 0, 1, ±N, N = 0,1, 2...) 
and corresponding two-component wave functions ^ n k 
are found from the Hamiltonian H (1) using the Landau 
gauge A = (0, Bx) and raising and lowering operators 
a + = Ibtt + /V2 and a — 1^77 /\/2, as in the work^. The 
magnetic length Ib and the c yclotron frequency u> c are 
defined as usual: Ib = y/h/eB, lo c = eB/m, e is the elec- 
tron charge. The basis consisting of the wave functions 
describing the states in the ordinary two-dimensional 



electron gas 



e v 4>Nk(x) is used, where k is the 



H — Hq + Hi + H2 



(1) 



parameter which labels degenerate states within one Lan- 
dau level in Landau gauge. 



Eo(t) = hu, E 1 @ = hu-t6+(—-£-)tooc (5) 
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J_0//7r+7r 
2m 71 V -7T7H 
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(3) 




ihky are the com- 



plex momentum operators, k is the wave vector measured 
from the center of the valley, £ is the valley index, £ = 1 
in the valley K, £ = —1 in the valley K' , 70 is the intra- 
layer A-B coupling parameter, 71 is the inter-layer A2-B1 
coupling parameter, m = r fi/2v 2 is the effective mass for 

bilayer graphene, v = -57^70, a is the lattice constant. 

The parameter 74 describes the next-nearest-neighbor 
transitions (A1-A2 and B1-B2 interlayer hopping), 71 = 
O.I70, 74 = 0.05 7o (see££). 

Ho is the basic term yielding a parabolic spectrum with 
the effective mass m, and its quasiparticles are chiral with 
the degree of chirality related to Berry phase 2n. 

Hi describes the layer asymmetry, leading to the open- 
ing of a gap ~ U in the spectrum. 

H 2 is due to the next-nearest-neighbor transitions and 
the difference between on-site energies of dimer and 



E±n(0 = ±hu c ^N{N - l)--^ + (i±-_)ft Wc (27V-l) 

2 70 271 

(6) 

where 5 — Uhuj c /ji. 



*ofc = (<?ta,0), *ife = Oifc,0), 



^nfcj = {a<nZ<t>Nk,b n £(f)N-2,k) 



(7) 



The coefficients a„j and b n £ are the eigenvector com- 
ponents. 



a n£ = 1/V 1 + D n f, b„f = D n f/Jl + D„ 



D 



E. 



g - jU/2 + £N5 - - ^)huj c (2N - 1) 
fiu e y/N(N-l) 



(8) 

Without any asymmetry in the zero approximation 
a±iv,s = l/V2,b ±Nti = ±l/V2. 

The asymmetry splits the zero-energy Landau level de- 
generate in valleys and orbital momenta into four levels. 
Energy levels for N > 2 are weakly split in valleys. 

Note that the spectrum of high-energy LLs is appli- 
cable for the fie lds and levels satisfying the condition 
huj cX /N(N - 1) < 71/4. For 71 = 0.39eV this inequality 
yields B < 50T for N = 2. For higher fields or higher 
levels the full four-band Hamiltonian has to be used to 
determine the exact LL spectru m 31 ! 32 . 
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The Zeeman splitting is omitted, and all levels are dou- 
bly degenerate in spin. Although in graphite the electron 
g-factor is not small (g = 2), a very light effective mass 
m w 0.054 in the bilayer determines a small ratio between 
the Zeeman energy and LL splitting ezjfy^c ~ 0.05 9 . 
Trigonal warping coming from 73 = 7A1-S2 *C 71 is not 
included. 



III. COULOMB INTERACTION AND 
MAGNETOEXCITATIONS 

The total Hamiltonian of the many-body system in the 
perpendicular magnetic field with the Coulomb interac- 
tion is 

H = E n £a^£ a a\£0 + H int (9) 

where a^ CT and a\^ a are the one-particle creation and an- 
nihilation operators; A = (n,k), n — 0,1, ±7V indicates 
the Landau level; k is the parameter which labels degen- 
erate states within one Landau level in Landau gauge; £ 
and a are valley and spin indexes. 

H in t = 2^2 V \3,Xi a ti^ a X3i'a' a ^e'r' a \l^ ( 10 ) 

The Coulomb interaction conserves spin and valley in- 
dexes; fcl + k2 = k3 + fc4. 

The matrix elements for Coulomb interaction are found 
using the two-component wave functions (7), in analogy 
with calculations for monolayer graphene^. 

= V^Y^'^'^^uMJu^-^ (11) 

(q) + (q) (12) 

■Wo) = (£r) 1/2e ^(^7f^) m "" L ^"(y)' ( 13 ) 

where V(q) = f|, q = (q x ,q y ),k4 = kl + q y , fc3 = k2 - 

%■ Jm,n(q) = ^,m(-q)( TO > n )'> L n~ n are Laggerre 
polinomials. 

In Eq. (10)the summation is over the ensemble 
nl,n2,n3, nA, fcl,fc2, both spins, both valleys and the 
wave vector q = (q x , q y ). 

In this work only the charge-density-excitations (tran- 
sitions without changing the electron spin) are studied, 
valley and spin indexes (£, a) are not changed. Corre- 
sponding operators for excitations (n,n') from the level 
n to the level n' with the momentum K are 

C'^W^E^W (14) 

k 

where A = (n, k), A' = (n' , k + K). It is assumed that the 
magnetic field is high which means that E c <C Tilo c , where 



E c is the typical Coulomb energy: E c — e 2 /els- The 
momentum of excitation is small: KIb <SC 1. The prob- 
lem is considered in the way analogous to that employed 
ir>22 for the two-dimensional electron gas with quadratic 
dispersion and i n 19 i 20 ' 25 for monolayer graphene systems. 
The time-dependent Hartree-Fock approximation is used. 
The Hartree-Fock approach assumes that there is a small 
parameter E c /AE nn t(£) <C 1, where AE nn /(£) is the 
transition energy without interaction. 

AE nn ,(0=E n ,(0-E n (0 (15) 

For monolayer graphene the ratio AEio/E c — 2.77 
(seei£) and it does not depend on the magnetic field 
strength due to linear dispersion. As stated in2£, the 
method works better for spin-flip excitations. For bilayer 
graphene E c = IOa/S, hu> c = 2.25 (see£) and the ratio 
hu} c /E c = 0.22S 1 / 2 for e = 5. For the first high-energy 
transition E12 ~ \/2hu} c , and therefore for B — 40T the 
ratio Ei2/E c ~ 2. We do not consider with this method 
the low-energy transitions between Landau levels and 
1 with energies close to zero. 

The excitation energy E n ^ n i.^ a consists of noninteract- 
ing and Coulomb parts: 

E n ,n>;ta = AE nn i (£) + E^ n , (16) 

Coulomb part E nn , is represented by the three terms: 
"excitonic" E^ n , part due to direct interaction of the 
electron at the level n' and the hole at the level n, 
exchange self-energy £ n f CT and S n '^ CT corrections to the 
one-electron Landau level energies and "depolarization" 
shift which is given in the random phase approximation 
(RPA). The RPA part is proportional to K, and is impor- 
tant for dispersion (Ref* 2 ^,). However, we are interested 
in the terms responsible for intervalley splitting which are 
constants independent of K. Therefore the RPA part can 
be omitted. This restriction enables to consider excita- 
tions with different (£, a) independently. 

En,n';£<7 = -^nfn' + ^n'£cr — (17) 

As for monolayer graphene, there is the problem of di- 
vergency of exchange self-energy E„ due to summation 
over all filled LLs. The spectrum of monolayer and bi- 
layer graphene described by the model Hamiltonian in 
unbounded both from above and below. This fact is 
physically artificial. In 2 ^ the cut-off value on energy or 
number of LL vas defined. Ir*i£ the semi-empirical way 
was used to treat the problem. The electron-electron 
interaction parameter is fitted for one type of transi- 
tion to experimental data. For the bilayer graphene the 
area of parabolic dispersion is less than required cut- 
off value, and it is necessary to consider the four-band 
Hamiltonian 2 ^—. In this work some general rules are 
presented allowing to conclude where the Coulomb part 
can be seen. 

The interlayer electron transitions from the top filled 
(fully or partially) to the next free (fully or partially) 
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Landau levels with energies nearly ui c are considered. 
The selection rules for these transitions are AN = 1. 
The case of filling-factor v — is considered. The Fermi 
level is equal to zero. This filling-factor means the ab- 
sence of charge density (the absence of free carriers or the 
equal amount of holes and electrons). Different possible 
ground states in magnetic field and different cyclotron 
transitions may correspond to this filling. 

A. The asymmetric bilayer without e-h asym- 
metry. 

Let U > 0. In this case we have filling- factor v = 4 
for the electrons in one valley and v = 4 for the holes 
in another valley. For the valley with £ = 1 there is the 
top filled LL with n = —2 and the transition (-2,1), and 
for the valley with £ = — 1 there are the top filled and 
1 LLs and the transition (1,2). Including spin there are 
two transitions of each type. The noninteracting part is 
the same for both types of transitions: 

AE_ 2A (l) = AE li2 {-l)=u c V2 + hu-5\ (18) 

Note that the electron-hole symmetry of one-particle 
Hamiltonian leads to the fact that (-2,1) and (1,2) tran- 
sitions are really the same and have the same energy. 
(-2,1) in electron representation is (1,2) in hole represen- 
tation. Taking into account spin degeneracy we have four 
transitions with equal energies. Since all types of asym- 
metry and electron-electron interaction are considered as 
small perturbations, the many-body corrections can be 
calculated with symmetric wave functions. 

For v = integer filling, labelled 0/, 
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dqy(g)J 22 (q)Ji,i(-q) 



2(2vr) 



dqV(q)(J 22 (q) + J 00 (q))Jll(q)) (19) 



(£ a - SOo/ = J dq^jjd Jll(q)| 2 + |Jlo(q) 
-/ dq 2W^ |J2l(q)|2 + ^ J20(q)|2)+ 




60 80 
Magnetic field (H) 



FIG. 1: The one-electron energy levels and cyclotron tran- 
sitions taking into account the next- nearest-neighbor transi- 
tions and layer asymmetry (solid for K valley and dash for K' 
valley) . 



and the other <j 2 is completely empty in both valleys. 
There are two excitations Qi 2 .^ ai and two excitations 
Q—2 if CT 2 with the same energies. 

B. The asymmetric bilayer with electron-hole 
asymmetry. 

In this case at the total filling-factor v = the filling 
of Landau levels depends of the magnetic field. 
If the magnetic field is not so high and 



1± _ J± 

7o 27; 



■)TlU! c — s < 



u 



(21) 



then in the valley with £ = — 1 levels and 1 arc filled and 
in the other valley with £ = 1 they are empty. There are 
cyclotron electron or hole transitions (1,2) from the top 
filled to the next empty level as for the case A (see Fig.l, 
left). The noninteracting parts of transition energies are 
presented by the following expressions: 



A 



A£ 1)2 (-l) = Huj c (V2 + - — ) 
7o 271 



5) (22) 



+ J ^^2- E l J 2,iY(q)Jo>- 2 q| (20) 

This value may depend on resolving the divergency prob- 
lem, but it is not zero. There is no Kohn's theorem^! 

For the excitations in the different valleys energy split- 
ting due to the layer asymmetry is absent. 

The same transitions may occur for the spin ferromag- 
netic state, where Zeeman splitting would be included: 
one spin component o\ of LLs and 1 is completely filled 



AE_ 2A (l) = nw c (V2-^ + ^) + hu-6) (23) 
70 271 2 

The Coulomb parts are the same, as it was discussed 
before. The difference in the energies between the valley 
transitions is due to the electron-hole asymmetry. 

A£7i, a (-1) - A£L 2il (l) = - -)hcj c (24) 

70 71 
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If the magnetic field is sufficiently high and 



( 7T- ) huj c 

7o 271 



x U 
S> 2 



(25) 



then in both valleys levels 1 are empty, but the level is 
filled for the K valley and empty for the K' valley (see 
Fig.l, right). There are only hole- type transitions (-2,1). 
The following noninteracting parts of transition energies 
are 

AE_ 2>1 (l) = nu, c (V2-^ + A.) + hu-S) (26) 
7o 271 2 



AE_ 2il (-l) = hcu c (V2- 



74 
7o 



A 

271 



+ t^)-t;(U-5) (27) 



Because of different filling of Landau levels in the val- 
leys the influence of Coulomb interaction differs in the 
part of the self-energies. 
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FIG. 2: The one-electron Landau energy levels and cyclotron 
transitions for symmetric and asymmetric (B) cases. For the 
symmetric case there are two transitions in each valley. 



~£ 2 Q-£io = /dq-^4(kio(q)| 2 -^2o(q)| 2 ) (28) 



Ei °4 S2o= iW! (29) 



B_ 2ll (l) - £- 2 ,i(-l) = (U-8) + ^c^f (30) 

The splitting between valley cyclotron transitions is 
due to Coulomb interaction and layer asymmetry 

C. The bilayer graphene in the full-symmetric 
case. 

In this case labelled h we have two half-filled zero- 
energy levels in both valleys: {uq — 1/2) and \{yi = 
1/2); this means v = 2 for the electrons in each valley and 
v = 2 for the holes in each valley. For each valley there 
are two transitions: (1,2) from half-filled to empty levels 
and (—2,1) from filled to half-empty levels. These transi- 
tions are connected by the Coulomb interaction V. Using 
the Hartree-Fock approximation for non-integer filling- 
factor o 35 i 36 two combined modes Qf a with the energies 
E s a are found. 



Qi 



1 

V2 



E a ,a=u e V2 + E^±-V 
2 



(31) 



These modes may be called symmetric and antisymmet- 
ric in analogy to modes in semiconductor bilayer. 



E. 



1 



-E 



(1,2) 



1 



ex,oi t ( s 2 (E 2 - 7^ -(S 2 - Si) J 



T> _ T/1,2 _ T/-24 _ 

K — v -2.1 — "l,2 — 



(2^) 2 



dq^(g)Ji 2 (q)J-2,i(-q) = 



w/ W(9)|Ju|a(q) = 5^*^ 



(33) 



(32) 



V" = 2.5^ and for B = 40T V" ~ 15meV. 

This splitting for combined electron-hole transitions 
from half-filled level is not specific to bilayer graphene. 
In monolayer graphene with valley asymmetry consid- 
ered ii*i£ the filling-factors were v = 2 for electrons in 
one valley and v = 2 for holes in another valley. In the 
simple case with symmetric valleys we have the half-filled 
zero-energy level in both valleys which means v = 1 for 
both electrons and holes in each valley. For monolayer 
graphene with half-filled zero-energy level there are re- 
lated transitions (0,1) for the electrons and (—1,0) for 
the holes, and the corresponding value of splitting for 
the combined modes V mg — \\f\E c is found using the 
wave functions from 1 ^. This value is nearly the same as 
for bilayer graphene (V mg ~ 2.5\^B), but for monolayer 
graphene it is possible to observe this splitting for lower 
experimentally used magnetic fields. 

D. The symmetric bilayer with electron-hole 
asymmetry. 

For the bilayer graphene with symmetric layers and 
therefore the symmetric valleys, but with the electron- 
hole asymmetry included there are half-filled Landau 
levels (z^o = 1/2) and empty 1 Landau levels (z/1 = 0) in 
both valleys. For each valley there are only hole-type 
(—2, 1) transitions. Coulomb corrections are equal for 
both valleys due to the same filling, and there is no split- 
ting in this case. 



G 



IV. SUMMARY AND DISCUSSION 

In conclusion, the cyclotron transitions for clean bi- 
layer graphene are studied. The electron-hole asym- 
metry, the small layer asymmetry and the influence of 
Coulomb interaction are taken into account. The charge- 
density excitations at small momenta are considered 
within the Hartree-Fock approximation. It is shown that 
the energies of cyclotron transitions in two valleys can 
be equal or can be split either due to the electron-hole 
asymmetry or due to Coulomb interaction. The splitting 
depends on applied magnetic and electric fields. At lower 
magnetic fields the energy splitting in the spectrum is 
due to electron- hole one-particle asymmetry while higher 
magnetic fields the energy splitting in the spectrum is 
due to Coulomb interaction. For the fully symmetric 
case with half-filled zero-energy levels (where electron- 
hole and layer asymmetries are not taken into account) 
the energy splitting proportional to Coulomb interaction 
is found both for bilayer and monolayer graphene. 

In the work 2 — the many-body corrections to cy- 
clotron resonance in monolayer and bilayer graphene 
due to Coulomb interaction were studied neglecting the 
electron-hole asymmetry. In the following work^S of the 
same author the weak electron-hole asymmetry was in- 
cluded, where the cyclotron resonances were considered 
using the four-band Hamiltonian. The four-band Hamil- 
tonian consideration is more precise, though the results 
concerning the noninteracting part coincide in the lim- 
its of accuracy. The employed SMA method to study 



many-body corrections due to Coulomb interaction (sin- 
gle mode approximation) for integer filling-factors is iden- 
tical to that used in the present work. Unfortunately it is 
impossible to compare directly the results of these work 
with ours in the aspect of many-body corrections because 
i n 29 ' 30 the different filling-factors (y = 4, 8, 12, 16) were 
considered. 

In the experimental work 7 - graphene bilayers embedded 
in the multilayer epitaxial graphene were studied. The 
splitting between electron-type and hole-type transitions 
is found in a relatively narrow range of B. The possible 
explanation invoking the electron-hole asymmetry yields 
the right values, but should be seen for wide range of 
magnetic fields. Moreover, in that work the other higher 
energy transitions corresponding to An = 1 were studied. 
In our notation the experimentally studied transitions are 
(-2,3) and (-3,2). To compare the experimental results 
with the theory it is necessary to study these transitions 
in the future. 

In this work there are a simultaneous consideration of 
both Coulomb interaction and electron-hole asymmetry, 
and suggestions concerning observation of splitting due 
to these factors at different regimes. 
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